Krylov Approximation of Linear ODEs with Polynomial Parameterization
نویسندگان
چکیده
منابع مشابه
Krylov Approximation of Linear ODEs with Polynomial Parameterization
We propose a new numerical method to solve linear ordinary differential equations of the type ∂u ∂t (t, ε) = A(ε)u(t, ε), where A : C → C is a matrix polynomial with large and sparse matrix coefficients. The algorithm computes an explicit parameterization of approximations of u(t, ε) such that approximations for many different values of ε and t can be obtained with a very small additional compu...
متن کاملPolynomial odes: examples, solutions, properties
ABSTRACT. Let P be a polynomial from R → R and D ∈ R. I will consider the properties of the class of ODEs Y ′ = P (y) ; Y (0) = D and their solutions. The solution space to these ODEs form a proper subspace of the analytic functions with ODEs. They have many interesting algebraic and topological properties. I will present efficient methods for generating the power series solutions to these poly...
متن کاملRelaxed Krylov subspace approximation
Recent computational and theoretical studies have shown that the matrix-vector product occurring at each step of a Krylov subspace method can be relaxed as the iterations proceed, i.e., it can be computed in a less exact manner, without degradation of the overall performance. In the present paper a general operator treatment of this phenomenon is provided and a new result further explaining its...
متن کاملGeometric integration of conservative polynomial ODEs ✩
We consider systems whose Hamiltonian is of the form H(q,p) = 2p + V (q), where the potential V is either cubic or quartic with no cubic terms. For most of these systems (in the measure sense) we give an explicit numerical integration method that preserves both phase space volume and the value of the Hamiltonian. This is exemplified in the Hénon–Heiles system. An application is to the hybrid Mo...
متن کاملNotes on Linear Odes
Remark 1.3. Suppose the characteristic polynomial is given by χA(λ) = Π m i=1(λ−λi) , where λi are distinct. Then the polynomials p`,j(t) have degree strictly smaller than νj (Exercise: Why?) If A is real, we want to find an expression for general real solutions. These are in particular the solutions one obtain when given real initial data. Before we proceed, let us just make a simple observati...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
ژورنال
عنوان ژورنال: SIAM Journal on Matrix Analysis and Applications
سال: 2016
ISSN: 0895-4798,1095-7162
DOI: 10.1137/15m1032831